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ABSTRACT 

We investigate the properties of r-mode and inertial mode of slowly rotating, non-isentropic, 
Newtonian stars, by taking account of the effects of the Coriolis force and the centrifugal force. 
The Coriolis force is the dominant restoring force for both r-mode and inertial mode, which are 
also called rotational mode in this paper. For the velocity field produced by the oscillation modes 
the r-mode has the dominant toroidal component over the spheroidal component, while the 
inertial mode has the comparable toroidal and spheroidal components. In non-isentropic stars 
the specific entropy of the fluid depends on the radial distance from the center, and the interior 
structure is in general divided into two kinds of layers of fluid stratification that is stable or 
unstable against convection. Because of the non-isentropic structure, low frequency oscillations 
of the star are affected by the buoyant force, which has no effects on oscillations of isentropic 
stars. In this paper we employ simple polytropic models with the polytropic index n = 1 as 
the background neutron star models for the modal analysis. For the non-isentropic models we 
consider only two cases, that is, the models with the stable fluid stratification in the whole 
interior and the models that are fully convective. For simplicity we call these two kinds of models 
"radiative" and "convective" models in this paper. For both cases, we assume the deviation 
of the models from isentropic structure is small. Examining the dissipation timescales due to 
the gravitational radiation and several viscous processes for the polytropic neutron star models, 
we find that the gravitational radiation driven instability of the nodeless r-modes associated 
with I' — \m\ remains strong even in the non-isentropic models, where /' and m are the indices 
of the spherical harmonic function representing the angular dependence of the eigenfunction. 
Calculating the rotational modes of the radiative models as functions of the angular rotation 
frequency fi, we find that the inertial modes are strongly modified by the buoyant force at small 
f2, where the buoyant force as a dominant restoring force becomes comparable with or stronger 
than the Coriolis force. Because of this property we obtain the mode sequences in which the 
inertial modes at large fl are identified as g-modes or the r-modes with I' = \m\ at small fi. We 
also note that as fl increases from 51 = the retrograde 5-modes become retrograde inertial 
modes, which are unstable against the gravitational radiation reaction. 



Subject headings: instabilities — stars: neutron — stars: oscillations — stars: rotation 



1. Introduction 

Since the recent discovery of the gravitational radiation driven instability of the r-modes by Andersson 
(1998) and Friedman & Morsink (1998), a large number of papers have been published on the properties 
of the rotational modes of rotating stars because of their possible importance in astrophysics (e.g., Kojima 
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1998, Andersson, Kokkotas & Schutz 1998, Kokkotas & Stergioulas 1998, Lindblom & Isper 1998, Beyer & 
Kokkotas 1999, Lindblom, Mendell & Owen 1999, Kojima & Hosonuma 1999, Lockitch & Friedman 1999, 
Yoshida & Lee 2000, Lockitch 1999, Yoshida et al 1999). Here, we have used the term of "rotational mode" 
to refer to both r-mode and inertial mode of rotating stars. As an important consequence of the instability, 
for example, Lindblom, Owen & Morsink (1998) argued that due to the instability, the maximum angular 
rotation velocity of hot neutron stars is strongly restricted. Owen et al. (1998) also suggested that the 
gravitational radiation emitted from hot young neutron stars due to the r-mode instability would be one of 
the potential sources for the gravitational wave detectors, e.g. LIGO. 

The dominant restoring force for the rotational modes is the Coriolis force and the characteristic 

frequency is comparable to the angular rotation frequency CI of the star (e.g., Greenspan 1964, Pcdlosky 
1987). The r-modes are characterized by the properties that the toroidal motion dominates the spheroidal 
motion in the velocity field, and that the oscillation frequency uj observed in the corotating frame of 
the star is given by a; = 2mSl/{l'{l' + 1)) in the Umit of f2 ^ 0, where m and I' are the indices of the 
spherical harmonic function Y™ representing the toroidal component of the velocity field of the mode (e.g., 
Papaloizou & Pringle 1978 , Provost, Berthomieu & Rocca 1981, Saio 1982). On the other hand, the 
inertial modes have the comparable toroidal and spheroidal components of the velocity field and therefore 
they arc not necessarily represented by a single spherical harmonic function (e.g. Lee, Strohmayer & Van 
Horn 1992, Yoshida & Lee 2000). Note that the inertial mode of this kind is also called "rotation mode" or 
"generalized r-mode" (Lindblom & Ipser 1998, Lockitch & Friedman 1999). 

Most of the recent studies on the r-mode instability driven by the gravitational radiation reaction 
are restricted to the case of neutron stars with isentropic structure, since the barotropic (one-parameter) 
equation of state, p = p{p) is usually assumed for both the interior structure and the small amplitude 
oscillation. In isentropic stars the specific entropy of the fluid in the interior is constant and the buoyant 
force has no effects on the oscillation. However, stars in general have non-isentropic interior structure, 
which may be divided into two kinds of layers with fluid stratification stable or unstable against convection. 
It is well know that the property of the oscillation modes of non-isentropic stars is substantially different 
from that of isentropic stars, since the buoyant force comes into play as a restoring force for low frequency 
oscillations like g-mode propagating in the layers with the stable stratification (see, e.g., Unno et al 1989). 
This is also the case for rotational modes, which have low oscillation frequencies comparable to the spin 
rate O of the star. For example, we note that the r-mode associated with 1/ = |m| having the nodeless 
eigenfunction in the radial direction is the only r-mode permitted in an isentropic star for a given m 
(e.g, Lockitch & Friedman 1999, Yoshida & Lee 2000), but in non-isentropic stars with the stable fluid 
stratification all the r-modes with I' > \m\ are possible for a given m, and there are in principle an infinite 
number of r-modes for given I' and m (e.g.. Provost et al 1981, Saio 1982, Lee & Saio 1997). We therefore 
consider it important to investigate the modal properties of the rotational modes in non-isentropic stars, 
particularly in order to answer the question whether the property of the r-mode instability for neutron star 
models with non-isentropic structure remains the same as that for isentropic neutron star models. 

In this paper, we study the r-mode and the inertial mode oscillations of slowly rotating, non-isentropic, 

Newtonian stars, where uniform rotation and adiabatic oscillation are assumed to simplify the mathematical 
treatments. The effects of the rotational deformation of the equilibrium structure on the eigenfunction 
and eigenfrequency are included in the analysis, using the formulation by Yoshida & Lee (2000) (see also 
Lee & Saio 1986, Lee 1993). We employ simple polytropic models with the polytropic index n = 1 as 
the background neutron star models. For the non-isentropic models we consider only two cases, that 
is, the models that have stable fluid stratification in the whole interior and those that are in convective 
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equilibrium. In this paper, following the conventions used in the field of stellar evolution and pulsation, we 
simply call these two kinds of models "radiative" and "convective" models, respectively, although no energy 
transport in the interior is considered. The plan of this paper is as follows. In §2, we briefly describe the 
basic equations employed in the linear modal analysis of slowly rotating stars. In §3, we show numerical 
results for the rotational modes of non-isentropic stars, and discuss the properties of the eigenfrequency 
and eigenfunction of the modes. In §4, we examine the stability of simple neutron star models against the 
r-modes, computing their growth or damping timescales due to the gravitational radiation reaction and 
some viscous damping processes. §5 is for conclusions. 

2. Method of Solutions 

The method of solutions used in this paper is exactly the same as that applied in Yoshida & Lee 
(2000). Here, therefore, we briefly describe the basic equations and the method of calculation. The details 
are given in Yoshida & Lee (2000). 

2.1. Equilibrium Configurations 

We consider oscillations of uniformly rotating Newtonian stars. Thus the angular rotation frequency of 
the star, Q is constant. We assume the polytropic relation for the equilibrium structure of the model, and 
the relation between the pressure p and the mass density p is given by 

p = Kp'+i, (1) 

where n and K are the polytropic index and the structure constant determined by giving the mass 
and the radius of the model, respectively. In this investigation, assuming slow rotation, we employ the 
Chandrasekhar-Milne expansion (see, e.g, Tassoul 1978) to obtain the equilibrium structure of rotating 
stars. In this technique the effects of the centrifugal force and the equilibrium deformation are treated as 
small perturbations to the non-rotating spherical symmetric structure of the star. The small expansion 
parameter due to rotation is chosen as Q = il{R^/GMy^'^, where R and M are the radius and the mass 
of the non-rotating star, respectively. The equilibrium configurations are constructed with accuracy up to 
the order of fl^. For simplicity, we consider a sequence of slowly rotating stars whose central density is the 
same as that of the non-rotating star. 

2.2. Perturbation Equations 

In the perturbation analysis, we introduce a parameter a that is constant on the distorted effective 
potential surface of the star. In practice, the parameter a is defined such that 

*(r,e) = *o(a), (2) 

where and ^ are the effective potentials for the non-rotating and slowly rotating models, respectively. 
Note that is equal to the gravitational potential and \1/ is the sum of the gravitational potential and the 
centrifugal potential. With this parameter a, the distorted equi-potential surface may be given by 



r = a{l + e{a,e,ip)} , 



(3) 
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where e is the function of the order of , and is obtained by giving the functional forms of ^ and ^Pq using 
the Chandrasekhar-Milne expansion (Yoshida & Lee 2000). Here we have used spherical polar coordinates 
(r, 6*, tp). However, we employ hereafter the parameter a as the radial coordinate, in stead of the polar radial 
coordinate r. 

The governing equations of nonradial oscillations of a rotating star are obtained by linearizing the 
basic equations used in fluid mechanics (e.g., Unno et al 1989). In the following, we use 5 and A to denote 
the Eulerian and the Lagrangian changes of physical quantities, respectively. Since the equilibrium state is 
assumed to be stationary and axisymmetric, the perturbations may be represented by a Fourier component 
proportional to e«('^t+™¥')^ where a is the frequency observed in an inertial frame and m is the azimuthal 
quantum number. The continuity equation may be linearized to be 

5p=-V.{pC), (4) 

where ^* is the Lagrangian displacement vector, and we have made use of 5v^ — i(a + mr2)^* = io;^* with uj 
being the oscillation frequency observed in the corotating frame of the star. Here denotes the covariant 
derivative derived from the metric . 

In this study, we consider only adiabatic oscillations. Thus, the perturbed equation of state is given by 

5p^Tp(^-^+CA^ , (5) 



where F is the adiabatic exponent defined as 



and Ai is the generalized Schwarzschild discriminant^ given by 

A, = -V,p~~ ^V,p, (7) 

and the positive and negative values of Ai indicate the fluid stratification unstable and stable against 
convection, respectively. 

The linearized Euler's equation is 

- uj'^i, + V, + (5$^ + + ^ V,;p^ A^ + 2iujne\/jLp, , (8) 

where $ denotes the gravitational potential, and is the rotational Killing vector, with which the 3- velocity 
of the equilibrium fluid of a rotating star is given as = f2(/3*. The perturbed Poisson equation is 

V,VM$ = 47rG'(5/9, (9) 

where G is the gravitational constant. 

Physically acceptable solutions of the linear differential equations are obtained by imposing boundary 
conditions at the inner and outer boundaries of the star. The inner boundary conditions are the regularity 



^ Since the Schwarzschild discriminant A; is proportional to V^s with s being the specific entropy, a model with ^4^=0 
throughout the interior is called isentropic model in this paper. Using the same barotropic (one-parameter) equation of state, 
P = P{p) foi" both the equilibrium structure and the perturbation leads to isentropic models. 
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condition of the perturbed quantities at the stellar center. The outer boundary conditions at the surface 
of the star are Ap/p — and the continuity of the perturbed gravitational potential at the surface to the 
solution of ViV*(5$ = that vanishes at the infinity. 

In order to solve the system of partial differential equations given above, we employ the series expansion 
in terms of spherical harmonics to represent the angular dependence of the perturbed quantities. The 
Lagrangian displacement vector, ^* is expanded in terms of the vector spherical harmonics as 

oo 

C =Y^aSi{a)Yr(B,^Y^\ (10) 

i>|m| 

{^'(«)^r,,+T.(a)^Fr,Je-', (11) 

^ oc 

^'' = -37? E {i^^(«)>^r.^-r,(a)sin^?r,r,4e-* (12) 
(Regge & Wheeler 1957, see also Thorne 1980). The perturbed scalar quantity such as 5$ is expressed as 

oo 

<5$= ^ <5$i(a)y™(0,^)e-*. (13) 

;>|m| 

Here, the expansion coefficients Tit are called "axial" (or "toroidal") and the other coefficients are called 
"polar" (or "spheroidal"). It is known that these two set of components are decoupled when il = 0. 

Substituting the perturbed quantities expanded in terms of the spherical harmonics into the linearized 
basic equations (|), (|), (|) and (|), we obtain an infinite system of coupled linear ordinary differential 
equations for the expansion coefficients. Notice that the system of differential equations employed in this 
paper is arranged to give eigenfunctions and eigenfrequencies consistent up to the order of for slow 
rotation. The explicit expressions for the governing equations are given in Appendix of Yoshida & Lee 
(2000). 

Since the equilibrium state of a rotating star is invariant under the parity transformation defined by 
9 TT — 9, the linear perturbations have definite parity for that transformation. In this paper, a set of 
modes whose scalar perturbations such as 6^ are symmetric with respect to the equator is called "even" 
modes, while the other set of modes whose scalar perturbations are antisymmetric with respect to the 
equator is called "odd" modes (see, e.g., Berthomieu et al. 1978). For positive integers j = 1, 2, 3, • ■ •, 
we have I = \m\ + 2j — 2 and I' — I + 1 for even modes, and I = |m| + 2j — 1 and /' = / — 1 for odd 
modes, where the symbols I and have been used to denote the spheroidal and toroidal components of the 
displacement vector ^\ respectively (see equations (p^)~(p^). Note that we do not use the term of the even 
(the odd) parity mode to refer to the polar (the axial) mode, although this terminology is used traditionally 
in rclativistic perturbation theory (e.g., Regge & Wheeler 1957). 

Oscillations of rotating stars are also divided into prograde and retrograde modes. When observed in 
the corotating frame of the star, wave patterns of the prograde modes are propagating in the same direction 
as rotation of the star, whereas those of the retrograde modes are propagating in the opposite direction to 
the rotation. For positive values of m, the prograde and retrograde modes have negative and positive w, 
respectively. 
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For numerical calculations, wc truncate the infinite set of linear ordinary differential equations to 
obtain a finite set by discarding all the expansion coefficients associated with I higher than Imax- This 
truncation determines the number of the expansion coefficients kept in the spherical harmonic expansion 
of each perturbed quantity. We denote this number as k^ax- Our basic equation becomes a system of 
4 X kmax-^^ order linear ordinary differential equations, which, together with the boundary conditions, are 
numerically solved as an eigenvalue problem of w = cr + mfJ by using a Henyey type relaxation method 
(e.g., Unno et al . 1989, see also Press et al. 1992). 

We cannot determine a priori kmax, the number of the expansion coefficients kept in the eigenfunctions. 
In this paper, we start with an appropriate trial value of A:„„x to compute an oscillation mode at a given 
value of Q. Then, we increase the value of fc„ax) until the eigenfrequency converges within an appropriate 
numerical error, where the mode identification is done by examining the form and the number of nodes of 
the dominant expansion coefficients. We consider that the mode thus obtained with the sufficient number 
of the expansion coefficients is the correct oscillation mode we want for the rotating model at Q. 

In order to check our numerical code, we calculate eigenfrequencies of the r-modes of the n — 3 and 
r = 5/3 polytropic model at = 0.1. In Table 1, our results are tabulated together with those calculated 
by Saio (1982) by using a semi-analytical perturbative approach. As shown by Table 1, our results are in 
good agreement with Saio's results. 



3. Rotational Modes of Slowly Rotating Non-isentropic Stars 

Because it is believed that isentropic structure is a good approximation for young and hot neutron 
stars, in this study we concentrate our attention to the case that the deviation of the non-isentropic models 
from isentropic structure is sufficiently small. In this case, we may assume that the adiabatic exponent F is 
given by 

In , ^ 

where 7 is a constant giving the deviation from the isentropic structure. By substituting equation ( p^ ) into 
equation (^, we obtain an explicit form of Ai as: 

= -7Vilnp. (15) 

In this study, we only consider oscillations of polytropic models with the index rt = 1. Positive and negative 
values of 7 are respectively for convective and radiative models. Isentropic models are given by 7 = 0. Since 
our formulation for stellar oscillations of rotating stars is accurate only up to the order of Ct^ , we shall apply 
our method to the oscillation modes in the range of f2 < 0.1. 

For convenience of later discussions, it may be useful to review the properties of r-modes and inertial 
modes of isentropic models with 7 = (see, e.g., Yoshida & Lee 2000). Because of the asymptotic behavior 
of the eigenfrequency of the rotational modes in the limit of ^ 0, it is convenient to introduce a 
non-dimensional frequency k as follows: 

Since the rotational deformation is of the order of il^, the effects on the frequency of the rotational modes 
appear as the terms of the order of Q^. Therefore, we may expand the dimensionless eigenfrequency k of 
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the rotational modes in powers of fl as follows: 

K = K0 + K2CP + 0{n'^) , (17) 

where 

= WTT) ^''^ 

for the r-modes associated with the indices I' and m of vector harmonics with the axial parity (see, e.g. 
Papaloizou & Pringle 1978). For the inertial modes, there is no analytic expression for ko in general. 

As suggested by Provost et al. (1981) and Saio (1982), the r-mode associated with /' = \m\ having the 
nodeless eigenfunction T;/ (r) is the only r-mode permitted for isentropic stars for a given m. In the limit of 
n ^ 0, we have k ^ kq = 2m/l'{l' + 1). Among the expansion coefficients of the eigenfunctions the toroidal 
component Tii~\^^ is dominating. For the inertial modes, the value of k in the limit of Q ^ is also finite 
but different from kq = 2m/l'{l' + 1). 

The inertial modes are characterized by the two angular quantum numbers and m, where the value 
of Iq is the maximum index I of spherical harmonics associated with the dominant expansion coefficients 
of the cigcnfimctions (see Lockitch & Friedman 1999, Yoshida & Lee 2000). Note that our definition of Zq 
is not the same as that of Lockitch & Friedman (1999), but the same as that of Lindblom & Ipser (1998) 
and Yoshida & Lee (2000). Note also that odd (even) parity modes have odd (even) values oi Iq — \m\, and 
that the odd parity mode having Iq — \m\ = 1 is the r-mode associated with V = \m\. The value of kq for 
the inertial modes depends not only on the indices m and Zqi but also on the equilibrium structure of stars 
(e.g., Yoshida & Lee 2000). For given values of Iq and to, the number of different inertial modes is equal to 
Iq — |to|(> 2). However, there seems to be no good quantum number to classify the different inertial modes 
associated with the same /o and m. We here introduce an ordering integer n to use the labellings such as 
in and cj„ for these inertial modes. Since the frequencies of the inertial modes with the same Zq and to are 
found in a sequence given by, in the case of m > 0, 

W-(io-|m|-i)/2 < • • • < w_i < < Wo < wi < • • • < a;(io_|„|_i)/2 (19) 
for odd parity modes, and 

W_(;„-|m|)/2 < • • • < W-1 < < Wl < • • • < UJ(l„-\m\)/2 (20) 

for even parity modes, we may define the ordering number n as an integer that satisfies the inequality 
\lo — |to| — l|/2 > |n| > for odd parity inertial modes, and |Zo — |to||/2 > |n| > for even parity inertial 
modes (see Lockitch & Friedman 1998, Yoshida & Lee 2000). We apply this classification to the inertial 
modes of both isentropic and non-isentropic stars in the following discussions. 



3.1. r-modes 

For the radiative models, all the r-modes associated with V > \m\ are possible for a given to, and there 
are in principle a countable infinite number of r-modes for given /' and to. The r-modes with the same /' 

and TO, may be classified by counting the number of nodes of the axial components T;/ in the radial direction. 
For them the values of K2 are different, but the values of kq are the same and equal to 2m/l'{V + 1). It 
is important to note that this classification of the r-modes is applicable only in the limit of O — > and 
K ^ Kq = 2m/l'{l' + 1), since the components associated with the other values of I' will come into play for 
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large O (see, e.g. Provost et al. 1981, Saio 1982). Let us introduce an integer k to denote the number of 
nodes of the toroidal eigenfunction T//, where the node at the center of the star is excluded from the coimt. 
Then, the r-modes are completely specified by the three quantum numbers m, V and k. In this paper, for 
given V and m, we employ the notation of r^. to denote the r-mode whose toroidal eigenfunction Tii{r) has 
k nodes in the radial direction. 

In Figures 1 to 3, the dimensionless eigenvalues n are given versus Cl for several low-radial-order 
r-modes with m = 1, 2, and 3 for the n — 1 polytropic model with 7 = — 10~^, where the r-modes with 
V = |m|, /' = |m| -t- 1 and I' = |m| -|- 2 are given in each figure. Modes with I' = \m\ + 1 are even parity 
modes, and those with I' ~ \m\ and V — \m\ + 2 are odd parity modes. In the figures, the eigenvalues k of 
the ZQ-modes with — \m\ + 2fc + 1, calculated at = 0.1 for the n = 1 isentropic model, are designated 
by the filled circles, where the integer k corresponds to the number in the subscript of the notation r^. 
Figures 1 to 3 clearly show that the qualitative behavior of the eigenfrequencies as functions of Cl is not 
strongly dependent on the azimuthal quantum number m. From these figures we also understand that 
the r-modes of the radiative model with 7 < are divided into three distinctive classes according to the 
asymptotic behavior of the eigenvalues at large Cl. The first class of the r-modes consists of the nodeless 
ro-modes associated with I' = \m\. The eigenvalue k of the r-modes of this class is well approximated by 
equations (17) and (18) at any values of Q examined in this paper. Note that the rp-mode with m = I' = 1 
is a peculiar mode in the sense that k is almost constant since K2 = and k ~ kq + 0{fl*) (see Saio 1982, 
Greenspan 1964). The second class consists of the r^-modes with V = \m\ and fc > 0. Equations ( p^ ) and 
( |l8| ) can give good values for the eigenvalue k of the r-modes of this class only when Cl is sufficiently small. 
As f2 increases, k tends to the io inertial modes with Iq — \m\ + 2fc + 1 of the isentropic model. Note that 
the iQ-modes are odd parity inertial modes (see equations (19) and (20)). We therefore consider that at 
large the r^-modes of this class become identical with the ZQ-modes with Iq = |rM| + 2k + 1, for which 
K is almost constant. The third class of the r-modes consists of the r^-modes with I' > \m\ and A: > 0. 
The eigenvalue k of the r-modes of this class is well represented by equations (17) and (18) only when is 
sufficiently small, and it becomes small monotonically as increases. The r-modes of this class may have 
no corresponding inertial modes of the isentropic model. 

As proved by Friedman & Schutz (1978), oscillation modes whose frequency a observed in an inertial 
frame satisfies the condition a{a + mil) = acu < are unstable to the gravitational radiation reaction. The 
condition for the instability may be rewritten as < k < m for retrograde modes with uj > and m > 0. 
Figures 1 to 3 show that the r-mode with V > 2 satisfies the inequality < k < to for the range of Vt 
examined in this paper. It is also found that the r^-modes with V = 1 and fc > also satisfy the inequality 
< K < TO for n > 0. The reason for the instability of the r^-modes with V = \ and fc > may be 
understood by noting that the polar components with the index I — 2 can be responsible for the emission 
of gravitational radiation. We may conclude that the rg-mode with V = \m\ = 1 is the only r-mode that is 
stable against the gravitational radiation reaction. 

In Figure 4, the eigenfrequencies k, of three r-modes with m = 2 are plotted as functions of I7I, where 
we have assumed H = 0.01. The three modes are respectively the rp-mode with /' = to = 2, the ri-mode 
with I' = m = 2, and the ro-mode with I' = 3 and to = 2. These three r-modes belong to the first, second, 
and third classes of the r-modes discussed in the previous paragraph. The solid lines and the dotted lines 
are used to denote the cases of negative and positive 7, respectively. Note that although we can obtain the 
three r-modes without any difficulty for the radiative models with negative 7, we cannot always calculate 
the r-modes of the second and the third classes for the convective models with positive 7. Let us first discuss 
the case of the radiative models. As shown by Figure 4, the eigenvalue k of the rp-mode with /' = |to| = 2 is 
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practically independent of I7I and has values approximately equal to kq = 2m/l'{l' + 1) = 2/3. However, k 
of the ri-modc with /' = |m| = 2 is dependent on |7|. At small I7I, the ri-mode has the eigenvalue close to 
Kq = 0.517, which corresponds to the inertial mode with Iq — \m\ = 3 of the isentropic model, but it tends 
to 2m/l'{l' + 1) = 2/3 as I7I increases. The ro-mode with I' = 3 and m = 2 also has the eigenfrequency k 
which depends on |7|. Although the value of k is close to kq = 2m/l'{l' + 1) = 1/3 at large I7I, it tends to 
zero as I7I 0. This behavior of k — » in the limit of |7| ^ is consistent with the fact that the r-modes 
associated with I' > \m\ do not exist in isentropic stars, for which 7 = 0. Let us next discuss the case of the 
convective models with 7 > 0. The eigenvalue k of the ro-mode with I' = \m\ = 2 is only weakly dependent 
on 7, as in the case of the radiative models. Note that the curves for the ro-mode with I' = \m\ = 2 for 
negative and positive 7 are almost indistinguishable in Figure 4. The ri-mode with I' = \m\ =2 for the 
convective models behaves differently than that for the radiative models as I7I increases. In fact, the value 
of K decreases rapidly with increasing 7 and the ri-mode cannot be found at large values of 7 > 0. No 
ro-mode with I' = 3 and m = 2 are found for the convective models. 

Let us discuss the properties of the cigcnfunctions of the r-modes for the radiative models. Since the 
amplitude of the displacement vector is much larger than that of the scalar perturbations, in the following 
discussions we may show only the displacement vector of the modes. In Figures 5 to 7, we show the 
expansion coefficients Si, Hi and zT;/ as functions of the fractional radius a/R for the r-modes with m = 2, 
where we have assumed 7 = —10""*, and the amplitude is normalized at the surface such that iTii — 1. 
Here, I = \m\ and I' = I + 1 for even modes, and I = \m\ + 1 and I' = I — 1 ior odd modes. The modes 
given in Figures 5 to 7 are respectively the ro-mode with I' = 2, the ri-mode with I' = 2, and the ro-mode 
with I' = 3. In each figure, the cigcnfmictions for the cases of = 10~^ and f2 = lO"""^ are shown in panels 

(a) and (6), respectively. The panels (a) of Figures 5 to 7 show that the toroidal component iT;/ of the 
r-modes is dominating the other expansion coefficients when O is small. However, as shown by the panels 

(b) , although zT;/ of the ro-mode with I' = 2 remains dominant even at f2 = 0.1, the toroidal components 
iTi' of the ri-modc with V = 2 and the rg-mode with Z' = 3 arc not necessarily dominating any more at 
n = 0.1. The properties of the eigenfunctions of the ri-mode with V = m = 2 dX large Ct may be regarded 
as those of the inertial modes. It is interesting to note that the displacement vector of the ro-mode with 
I' = 3 has large amplitude in the deep interior and almost vanishing amplitude near the stellar surface at 
large Vl. This suggests that because of the low frequencies u) in the corotating frame the buoyant force is as 
important as the Coriolis force to determine the mode character at large Q. 



3.2. Inertial modes 

For the radiative models with the polytropic index n = 1 and 7 = — 10~*, the eigenfrequencies 

uj = ijj / {GM / R^Y/"^ of low frequency oscillations are shown as functions of f2 for m = 1, 2, and 3 in Figures 
8 to 10, where even and odd parity modes arc displayed in panels (a) and (b), respectively. The straight 
dashed lines drawn in the figures are given by a) = and u) = mfl (see below). In the figures, the inertial 
modes of the n = 1 isentropic model computed at = 0.1 are designated by the filled circles. We note that 
the oscillation modes in a mode sequence are regarded as inertial modes when 51 is large, and as internal 
gravity (g-) modes when Cl is small. If we follow one of the continuous mode sequences from large Cl to 
small O, we may physically understand the behavior of the low frequency oscillation in the radiative models 
as a function of Cl. At large f2, the Coriolis force is dominating as the restoring force and the oscillation 
mode is identified as an inertial mode. As O decreases, however, the buoyant force becomes dominating the 
Coriolis force and the mode becomes identical with a g-mode in the limit of 11 ^ 0. Because of the effects 



- 10 - 



of the buoyant force in the radiative model, the frequency w of the inertial modes deviates from that of the 
inertial modes of the isentropic model as fl decreases. 

The modes within the two straight dashed lines depicted in Figures 8, 9 and 10 arc unstable against 
the gravitational radiation driven instability in the sense that they satisfy the frequency condition 
a{a + mfl) < (Friedman & Schutz 1978). We can see from Figure 8 that the ifc-modes with k ^ 
and m = 1 arc stable against the instability for all f2. Figures 9 and 10 show that as f2 increases all the 
retrograde (7-modes with |m| > 2, which arc stable against the instability when il ~ 0, become retrograde 
inertial modes, which arc unstable against the instability at large fl. 

There are several cases where we find it practically impossible to obtain completely continuous mode 
sequences between the g(-modes at small f2 and the inertial modes at large Q. In such sequences, there 
exist a domain of f2, in which there appear a lot of modes having a large number of nodes of the dominant 
expansion coefficients, and it is difficult to identify the same oscillation modes as those computed in the 
previous steps in 0,. In Figures 8 to 10, the domains of Q. with this difficulty are indicated by the dotted 
lines in the mode sequences. Unfortunately, we cannot remove the difficulty by increasing kmax- One of 
the reasons for the difficulty may be attributed to the method of calculation we employ in this paper. As 
discussed in Unno et al (1989), since we have truncated the infinite system of differential equations to obtain 
a finite system of equations, the numerical procedure used here inevitably contains a process of inverting 
matrices which are singular at some value of the ratio 2f2/a; 1, depending on the oscillation mode we are 
interested in. In this paper, wc assume that there always exist a continuous mode sequence between the 
(/-mode and inertial mode as a function of f2, even if the mode sequence contains a domain of fl in which 
the modes we are interested in cannot be obtained numerically. 



In order to find the connection rule between the g-mocles at small and the inertial modes at large 
ri, we tabulate first the frequencies of the 5-modes at H = for the radiative model with the index n = 1 
and 7 = — 10~^ in Table 2. To obtain a clear understanding of the global correspondence between the g-, 
r-, and inertial modes, with the help of Table 2 and Figures 8 to 10, we give schematic diagrams for the 
connection rules between the several low-radial-order modes as follows: 

For odd parity modes we have 



3.3. Connection rules between the (7-modes, r-modes and inertial modes 
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\m\ 
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rn 



+ 5) 
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+ 5) 
+ 5) 



(21) 
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(22) 



and for even parity modes we have 
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(23) 



91 = \m\) 
g-i{l = |m|) 



where the two corresponding modes occupy the same position in the two matrices connected with an arrow. 

Prom the diagrams, we may find the connection rule between the pife-modes with / and m and the 
inertial modes i±j with Iq and m: 

5±fc(? = P)^i±,(Zo = P+2fc), (25) 

where P = \m\ + 2j — 2 for even modes and P = \m\ + 2j — 1 for odd modes, and j and k are positive 
integers. The symbol g±k denotes the g-mode with k radial nodes of the expansion coefiicient Si, having 
positive and negative frequencies oj observed in the corotating frame of the star when Cl^O. 

The connection rule between the r-modes associated with /' = |m| and the inertial modes io with Iq 
and m may be given by 

r,_i(/' = |m|) ^ io{lo = |m| + 2j - 1) , (26) 

where j is a positive integer. Note that the io-mode with lo = |m| + 1 is exactly the same as the ro-modes 
with I' = \m\. To establish the connection rules given above, we have neglected the effects of avoided 
crossings between the modes associated with different values of I and/or k, since the mode property is 
carried over at the crossings. 
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(24) 
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4. Dissipation Timescales for r-modes of Non-isentropic Stars 

For the rotational modes of non-isentropic models, we have found that all the r-modes except the 
TQ-mode with V — \m\ — 1 are unstable to the gravitational radiation reaction in the sense that the 
frequency a satisfies the condition a{a + mfl) < 0. Let us determine the stability of the r-modes of the 
non-isentropic models against the gravitational radiation driven instability, taking account of the dissipative 
processes due to the shear and bulk viscosity, where we will follow the method of calculation almost the 
same as that employed in the previous stability analyses of the ro-mode with I' ~ \m\ and the inertial modes 
for isentropic models (Lindblom et al 1998, Owen et al 1998, Andersson et al 1998, Kokkotas & Stergioulas 
1998, Lindblom et al 1999, Lockitch & Friedman 1999, Yoshida & Lee 2000). 

The damping timescale of the rg-mode associated with V — \m\ for sufficiently small Cl may be 
estimated by (Lindblom et al 1998) 

1 1 /io^K\" 1 / T V / n 







f T \ 









+ 



1=2 ^■^'^ y^^pj \ttgp 



where p is the average density of the star. The dissipative timescales f are calculated for each oscillation 
mode in the limit of f2 — > so that they are independent of the rotation frequency f2. Here, the first, 
second, third and fourth terms in the right-hand side of equation ( |27| ) are contributions from the shear 
viscosity, the bulk viscosity, the current multipole radiation and the mass multipole radiation, respectively. 
The expression (^7|) of the damping timescales has been derived on the assumptions that the frequency of 
the modes is well represented by a linear function of f2 and the toroidal component of the displacement 
vector is dominating. Although the assumptions are reasonable for the ro-modes with I' — \m\ (see Yoshida 
and Lee 2000), they are not necessarily correct for the r^-modes with I' ~ \m\ and k > 0, for which we use 

1 1 flO^K\^ 1 f T \^ firCp 



/lO^ifV 




( T \ 


[ T ) ~ 







fs\Tj fB yio^K J V 

l+l oo / r^2 \ '+2 



T.j-[^J +E?-b^ ' (28) 



"^j,; V'G'p/ ^fDj\TTGpJ 

because these r-modes become inertial modes at large fi. The difference between these two expressions is 
found in the second terms in the right-hand side of the two equations. When Q is small, equation ( p8| ) 
gives the correct expression of for inertial modes of isentropic stars (See, Lockitch & Friedman 1999, 
Yoshida & Lee 2000). To determine the damping timescales f in ( p8| ) in the limit of f2 — *■ 0, we assume 
that the dissipative timescale f~^ depends on fl as afl^ + b for each dissipative process, where the two 
unknown constants a and b are determined by calculating f^^ and df~^ /dfl at fl = 0.1. The f~^ in the 
limit of ^ is given by f~^ — b. We note that equation ( p8|) is not necessarily a correct approximation 



for non-isentropic stars. We use the expression (28) simply because we want to give rough estimates of the 



dissipative timescales, in order to compare the timescales with those of the isentropic case. 

In Table 3, we tabulate the dissipation timescales f in the unit of second for the various dissipative 
processes for the rfc-modes with I' = \m\, where the radius and the mass of the n = 1 polytropic neutron 
star model at r2 = are chosen to be i? = 12.57km and M = IAMq, respectively. Note that in Table 
3 we have not included the r^j-modes with V > \m\, since the frequency lu becomes vanishingly small at 
large f2 (see Figures 1 to 3) and hence the magnitude of the gravitational radiation driven instability, which 
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is proportional to w, becomes quite small. For the ro-modcs with I' = |m|, in order to see the effects of 
the buoyancy force on the stability, the results for the four different values of 7 arc tabulated. As shown 
by Table 3, the gravitational radiation driven instability for the rg-modes with I' = \m\ is not affected by 
introducing the deviation from isentropic structure and remains strong even for the non-isentropic models. 
As for the r^-modcs with /' = |m| and fc > 0, it is found that the dissipativc timescalcs are comparable to 
those of the corresponding inertial modes (compare Table 4 of Yoshida & Lee 2000). Thus, we may conclude 
that the r-mode instability driven by the gravitational radiation reaction is dominated by the ro-mode with 
I' = \m\ = 2 both in isentropic and in non-isentropic stars. 

5. Conclusions 

In this paper, we have numerically investigated the properties of rotational modes for slowly rotating, 
non-isentropic, Newtonian stars, where the effects of the centrifugal force are taken into consideration to 
calculate the modes with accuracy up to the order of Estimating the dissipation timescales due to the 
gravitational radiation and the shear and bulk viscosity for simple polytropic neutron star models with the 
index n= 1, we have shown that the gravitational radiation driven instability of the ro-mode with /' = |m| 
is hardly affected by introducing the deviation from the isentropic structure, and that the instability of the 
ro-modes with I' = \m\ = 2 remains strongest even for the non-isentropic models. Note that the ro-modes 
with /' = |m| arc the only r-modes possible in isentropic stars and their modal property docs not very much 
depend on f2 and the deviation from isentropic structure. We thus conclude that the result of the stability 
analysis of the ro-modes with I' = \m\ for isentropic stars remains valid even for non-isentropic stars. 

We have numerically found that for the radiative models the inertial modes at large O become identical 
with the rjt-modes with I' = \m\ and fc > and g-modes in the limit of O ^ because of the effect of 

the buoyant force, where the integer k denotes the number of node of the eigenfunction with dominant 
amplitude. This behavior of the (7-modes or the inertial modes as functions of fl has been suggested by, for 
example, Saio (1999) and Lockitch & Friedman (1999). We find the connection rules between the rfe-modes 
with /' = |m| and fc > and modes at small 0. and the inertial modes at large Cl. We also suggest that 
the rfc-modes with > \m\ and fe > 0, whose frequency in the corotating frame becomes vanishingly small 
as n increases, have at large 51 no corresponding inertial modes of the isentropic models. 

Good progress has been made in the understanding of the properties of the r-mode instability due 
to the gravitational radiation since its discovery. Recently, for example, Yoshida et al (1999) have shown 
that the results of the r-mode instability obtained by using the slow rotation approximation, as employed 
in this paper, are consistent with the results obtained by calculating directly the r-mode oscillations of 
rotating stars without the approximation. Kojima (1998) and Kojima & Hosonuma (1999) have studied 
the time evolutional behavior of r-modes by using the Laplace transformation technique for slowly 
rotating relativistic models. Lockitch (1999) also obtained r-modes and rotation modes for slowly rotating, 
incompressible, relativistic stars. However, we believe it necessary to show that the gravitational radiation 
driven instability of r-modes remains sufficiently robust in more realistic and unrestricted situations, 
taking account of the effects of differential rotation, magnetic field and general relativity, for example. To 
investigate the r-mode instability in more general situations is one of our future studies. 

We would like to thank Prof. H. Saio for useful comments and discussions. S.Y. was supported by 
Research Fellowship of the Japan Society for the Promotion of Science for Young Scientists. 
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Fig. 1. — Dimensionless frequencies k = u/^ for the low-radial-ordcr r-modcs with m = 1 are plotted as 
functions of the dimensionless angular rotation frequency Cl = {GM/ R^)^^^ , for the n = 1 polytropic 
model with 7 = —10^"'. The solid, the dotted and the dashed curves are for the ro, ri, and r2-modcs, 
respectively. The attached label I' denotes the angular quantum number I' associated with the spherical 
harmonic function representing the toroidal component of the displacement vector. Note that the eigenvalue 
of the r-modes is given by kq = 2m/{l'{l' + 1)) in the limit of O ^ 0. The filled circles at O = 0.1 denote 
the eigenvalues k of the inertial modes calculated for the isentropic model with the index n = 1 and 7 = 
a.tn = 0.1. 




Fig. 2. — Same as Figure 1, but for m = 2. 
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Fig. 3. — Same as Figure 2, but for m = 3. 
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Fig. 4. — Dimcnsionlcss frequencies k for the low-radial-order r-modes with m = 2 are plotted as functions 
of I7I for the n = 1 polytropic model at = 10~^. The attached label I' denotes the index of the spherical 
harmonic function for the toroidal component of the displacement vector. The number of nodes of iT;/ in 
the radial direction is indicated by the subscript in the notation rfe attached to the mode sequences. The 
solid line and the dotted line are given for negative and positive values of 7, respectively. 
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Fig. 5. — The first expansion coefficients Si, Hi, iTii for the tq mode with I' = m = 2 are plotted against 
a/R for the polytropic model with the index n = 1 and 7 = — 10~*. Panels (a) and (6) are for the cases 
of fi = 10^^ and 0, ~ 10^^, respectively. The cigenfonctions are normalized so that T|'(a = R) = 1. The 
attached labels Si, Hi, and iTc denote the expansion coefHcients of the displacement vector. 
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Fig. 7. — Same as Figure 5, but for the ro-mode with /' = 3. 
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Fig. 8. — Dimensionless frequencies w = uj / {GM / R^Y^'^ for the low-radial-order g- and inertial modes with 
m = 1 are plotted as functions of f2 = Q/{GM/R^y/^, for the n = 1 polytropic model with 7 = — 10~^. 

The even and the odd parity modes arc shown in panels (a) and (fe), respectively. The attached labels 
g±k{l) < — > i±j{lo — Iml) indicates the connection rule between the (;-modes and inertial modes. The filled 
circles at n = 0.1 denote the eigenfrequencies ui of the inertial modes for the n = 1 isentropic model at 
0, = 0.1. The part indicated by the dotted line in the mode sequence is the domain of Cl, in which the modes 
belonging to the sequence cannot be calculated. The two dashed lines are given by a) = and u) = mQ. 
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Fig. 9. — Same as Figure 8, but for m = 2. 
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Table 1. (w/fi — Ko)/ri^ for r-modes of the n = 3 
polytrope with F = 5/3 



m 


I' 


Mode 


Present 


Saio 1982 


1 


1 


^•0 


0.000 


0.000 






n 


-0.040 


-0.040 






r-2 


-0.115 


-0.114 




2 


r-o 


-0.100 


-0.099 






n 


-0.191 


-0.190 






r-2 


-0.315 


-0.312 




3 


r-o 


-0.037 


-0.037 






n 


-0.063 


-0.061 






r2 


-0.097 


-0.096 


2 


2 


ro 


0.121 


0.120 






n 


0.024 


0.024 






r2 


-0.001 


-0.001 




3 


ro 


-0.037 


-0.037 






n 


-0.071 


-0.070 






r2 


-0.115 


-0.114 


3 


3 


ro 


0.174 


0.173 






ri 


0.064 


0.063 






r2 


0.028 


0.027 



'^KQ is defined as ko = 2m/(/'(i' + 1)). The eigenfrequencies are 
calculated at Q = 0.1. 



Table 2. The eigenfrequency uj for low-radial-order 
5-modes of the n = 1 polytrope with 7 = — 10~^ at 
= 



I 


91 


92 


1 


1.20 X 10-2 


7.69 X 10-3 


2 


1.68 X 10-2 


1.15 X 10-2 


3 


2.01 X 10-2 


1.44 X 10-2 


4 


2.26 X 10-2 


1.68 X 10-2 


5 


2.47 X 10-2 


1.88 X 10-2 


6 


2.65 X 10-2 


2.05 X 10-2 
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Table 3. Dissipative timescales f of r- modes'' for the polytropic neutron star models with the index n = 1. 



m mode 7 kq ^fl(s) fs{s) •rx|m|(s)'' rD,\m\+i{s) T'j,|m|+2(s) 



1 


ri 







0.691 


5.89 X 10^ 


9.23 X 


107 




-2.46 


X 


10^ 


-1.27 X 


108 






-10" 


-4 


1.000 


6.08 X 10^ 


9.04 X 


10^ 




-2.31 


X 


10^ 


-1.58 X 


108 


2 


ro 


10- 


-4 


0.667 


2.03 X 10" 


2.50 X 


108 


-3.31 X 10" 


-3.49 


X 


102 















0.667 


2.03 X 10" 


2.50 X 


10* 


-3.31 X 10" 


-3.49 


X 


10^ 










-10" 


-4 


0.667 


2.03 X 10" 


2.50 X 


108 


-3.31 X 10" 


-3.49 


X 


10^ 










-10- 


-2 


0.667 


2.10 X 10" 


2.54 X 


108 


-3.31 X 10" 


-3.49 


X 


102 








ri 







0.517 


6.47 X 10^ 


6.18 X 


10^ 


-1.31 X 108 


-8.39 


X 


10* 


-1.88 X 


10'^ 






-10- 


-4 


0.667 


6.79 X lO'' 


6.01 X 


107 


2.18 X 107 


-7.05 


X 


10* 


-2.06 X 


10^ 


3 


ro 


10- 


-4 


0.500 


6.64 X 10^° 


1.43 X 


108 


-3.17 X 10^ 


-1.88 


X 


10^ 















0.500 


6.63 X 10^° 


1.43 X 


108 


-3.17 X IQi 


-1.88 


X 


103 










-10- 


-4 


0.500 


6.64 X 10^° 


1.43 X 


108 


-3.17 X 10^ 


-1.88 


X 


10^ 










-10- 


-2 


0.500 


6.97 X IQi" 


1.46 X 


108 


-3.17 X 10^ 


-1.87 


X 


10^ 








ri 







0.413 


6.97 X 10^ 


4.78 X 


107 


-1.86 X 10^" 


-5.30 


X 


10^ 


-4.07 X 


10"^ 






-10- 


-4 


0.500 


7.42 X 10^ 


4.63 X 


107 


2.39 X 10^ 


-4.37 


X 


10^ 


-4.48 X 


10"^ 



^We present dissipative timescales only for those that are unstable to gravitational radiation reaction. 

''We present dissipative timescales f; due to the gravitational radiation reaction only for the dominant 
multipole moments. 

"^For these dissipative timescales, we obtain positive values, which should be negative by definition. 
The reason for the inconsistent result may be because the extrapolation formula (p8|) for t~^, which is 
obtained at Q = 0, does not necessarily work well for arbitrary values of Q. For instance, we find that 
''^J,\m\/^D,\m\ + i ^ tO^ foi' rn = 2 and 'Tj^imi /t^j ~ 10^ for m = 3, where the dissipative timescales are 

estimated at f2 = 0.1. This may suggest that the mass multipole radiation dominates the current one for the 
ri-modes for large fi. Thus, our conclusions on the stability of the ri-modes extrapolated for large f2 are not 
altered at all. 



